Consistent picture for the electronic structure around a vortex core in iron-based superconductors 
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Based on a two-orbital model and taking into account the presence of the impurity, we studied theoretically 
the electronic structure in the vortex core of the iron-Pnictide superconducting materials. The vortex is pinned 
when the impurity is close to the vortex core. The bound states shows up for the unpinned vortex and are 
wiped out by a impurity. Our results are in good agreement with recent experiments and present a consistent 
explanation for the different electronic structure of vortex core revealed by experiments on different materials. 
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The iron-based superconducting (SC) materials have at- 
tracted much attention since their discovery JTl. Similar to the 
cuprate materials, the superconductivity in iron-pnictides can 
be achieved by doping either holes or electrons into the parent 
compound, among which the compound BaFe2As2, named as 
"122" system, is one of the most studied iron-based systems. 
The superconductivity can be realized by substituting the Bar- 
ium ions by Kalium ions or iron ions by the cobalt ions, cor- 
responding to the hole-doped materials Bai_ I K^FeiAs2 and 
electron-doped materials BaFe2- A Co A As2, respectively. 

The scanning tunneling microscope (STM) experiment has 
been a powerful tool to investigate the electronic structure of 
SC materials. Moreover, it can reveal the vortex core struc- 
ture and possible vortex state in presence of the magnetic 
field. Recently, the STM experiments were performed on 
both BaFe2- x Co A .As2 |2j] and Bai_^K x Fe2 As2 [3] compounds. 
However, the measurements on these two doped materials re- 
vealed significantly different results. For the hole-doped com- 
pound Bai_ v K A Fe2As2, the vortex forms a regular distorted 
triangle lattice. The strong bound state at the vortex core cen- 
ter was observed clearly with the energy peaked at the nega- 
tive energy. The intensity decreases and the peak splits with a 
dominant weight at negative energy when away from the cen- 
ter and finally evolves to the gap edges. However, for the 
compound BaFe2- A Co A As2, different results were reported, 
namely, the vortex forms disordered lattice and no bound state 
was observed at all. 

The contrast experimental observations in Bai_ A K v Fe2As2 
and BaFe2- A Co x As2 are intriguing. Theoretically, the exis- 
tence of bound states is robust and can be reproduced based 
on various effective models for iron-pnictides, while the po- 
sition of the resonance peak is model dependent 10-0] • Es- 
pecially, the negative bias peak is obtained ((J 0] based on a 
two-orbital model proposed in Ref. ]8|], a three-orbital model 
in Ref. [9], and a five-orbital one in Ref. [10]. Although there 
has been no theoretical calculation to explain the absence of 
the bound states in BaFe2- A Co A As2, a possible explanation is 
proposed in Refs. [0, @], i.e., the strong spin-density-wave 
(SDW) induced by the magnetic field inside the core would 
suppress the bound state. While the in-gap peak structure is 
still visible in the numerical results even in presence of the 
strong SDW order H. We note that for Cobalt doped sam- 



ples, the static SDW order disappears at the doping around 
6 = 0.065 ITT1 [l2tl and the STM experiment was performed 
at the doping density 6 = 0.1 ||2|], thus there is no SDW order 
in the bulk. The SDW fluctuation inside the vortex could be 
induced by the magnetic field and we expect it may suppress 
the bound state more or less while it is not strong enough to 
wipe out the bound states completely. On the other hand, this 
issue was also discussed qualitatively in Ref. [3] when com- 
paring their experimental results with Ref. [2]. It was pro- 
posed that the in-plane Cobalt ions or multi-band nature of 
the iron pnictides may account for such discrepancy. Since 
the Bai_ A K A Fe2As2 and BaFe2- A Co A As2 share the same par- 
ent compound, they should have similar multi-band structures. 
Thus we expect that the Cobalt ions, which may serve not only 
as dopants but also impurities, should be responsible for this 
discrepant behavior. In fact, the impurity effect and its role in 
the uniform SC state in iron-pnictides have been studied in- 
tensively dH [l3l - tl9ll and it was proposed that it can be used 
to probe the pairing symmetry Igl (TH HH. In presence of the 
magnetic field, the interplay between the impurity and vortex 
is of interest and studied in cuprates ll2"bl I2H1 . while the rele- 
vant study about the iron-pnictides is still awaited. Therefore, 
it is timely and interesting to address the impurity effect in 
presence of the field and present a consistent picture for the 
STM results for Co-doped and K-doped materials. 

In this paper, we study the single impurity effect based on 
a two-orbital model [8] which has given a sound explanation 
for the vortex bound state in Bai_ A K v Fe2As2 [6]. The model 
is given by 
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where i = (i x , i y ), j = (J x , j y ) are the site indices, yu, v = 1,2 are 
the orbital indices, and /i is the chemical potential. n- vcr is the 
density operator at site i and orbital /j, with spin cr. The quan- 
tity U' is taken to be U-2Jh- V s is the impurity potential. The 
impurities are located at the site i m . Generally speaking, for 
the magnetic impurity, the impurity scattering would include 
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both the potential scattering and the magnetic part. Based on 
first-principles electronic-structure calculations, the potential 
scattering potential is much stronger than the magnetic scat- 
tering one Il22ll . In the present work, the Cobalt irons are taken 
as a strong potential scattering center, while the weaker mag- 
netic part is not considered. 

The hopping constants are written as = 
%jv,o exp[/(;r/0o)] J*' A(r)dr in presence of the magnetic 
field, with <D>o is the superconducting flux quantum and 
A = (-By, 0,0) is the potential in the Laudau gauge. 
Following Ref. |8y, we use 
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The mean-field Hamiltonian can be diagonalized by solving 
the Bogoliubov-de Gennes (BdG) equations, 
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Here V; M j v is the pairing strength and f(x) is the Fermi-Dirac 
distribution function. 

The LDOS is calculated according to 

pAfJ) = YjKofW* -oi) + y^\ 2 S(E„ + o>)] , (10) 



where the delta function 8(x) is taken as F/n(x 2 +T 2 ), with the 
quasiparticle damping F = 0.01. 

In the numerical calculations, following previous stud- 
ies JHHH, we set h- 4 = 1,0.4,-2,0.04. The on-site 
Coulombic interaction U and Hund's coupling Jh are taken 
as 3.4 and 1.3, respectively. The pairing is chosen as next- 
nearest-neighbor intra-orbital pairing with the pairing strength 
V = 1.2, which will reproduce the s x 2 v 2-wave pairing symme- 
try, being consistent the experimental ll26l l27ll and theoretical 
proposal [28-31] for the gap symmetry in iron pnictides. The 
doping density is fixed at the electron-doped 6 — 0.1 accord- 
ing to the experimental value J2J]. A 40 x 20 square lattice 
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FIG. 1 : (Color online) The intensity plots of the gap magnitude in 
presence of one impurity at the site (10, 14), (13, 14), (16, 14) and 
(20, 10), respectively. The impurity locations are marked by the 
black spots. 
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FIG. 2: (Color online) The intensity plot of the magnetic order in 
presence of one impurity at the site (10, 14). 



with the periodic boundary conditions is used and a 10 x 20 
super-cell is taken to calculate the LDOS. 

The magnitudes of the SC gap with a strong impurity V s = 
100 at different locations are plotted in Fig. 1. Without the 
impurity, the vortex cores would appear at the site (10, 14) 
and (30, 6) (not presented here). When the impurity is put at 
the core center or near the vortex core, the vortex would be 
pinned by the impurity, as seen in Figs. 1(a) and 1(b). Here 
the pinned effect in the present two-orbital model is similar 
to previous results in cuprates 1I20I l2lll . The critical distances 
r c between the impurity and the core center (in a clean sys- 
tem) can be defined below which the vortex core center will 
be dragged to the impurity site. For the present parameters 
r c is about four lattice constant. It is insensitive to the direc- 
tion of the impurity but sensitive to the impurity potentials 
and electron interaction, which is not concerned with in the 
present work. As the distance is larger than the critical one, 
as seen in Fig. 1(c), the vortex and the impurity will be sepa- 
rated. Since the distance is just slightly larger than the critical 
one, for this case the vortex is still dragged to be close to the 
impurity site. The above pinned and dragged effects present a 
reasonable explanation why the vortex lattice is disordered in 
Co-doped materials |2J. As the distance increases further, as 
seen in Fig. 1(d), the impurity has no effect on the vortex. 

The intensity plot of the magnetic order (defined as 
1 1 /2 Zu(«i/i? - ijjuj.)!) is plotted in Fig. 2. As seen, the mag- 
netism is induced by the field and reaches the maxima value 
near the unpinned vortex core center. Away from the vortex, 
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FIG. 3: (Color online) The averaged low energy p(co)da) with 
e = 0.05) LDOS map with one impurity at (10, 14). 



the magnetism disappears gradually. These features are qual- 
itatively the same with previous studies about field induced 
SDW order without considering the impurity effect Jj-lal- For 
the case of pinned vortex, the magnetic order will be sup- 
pressed significantly near the impurity site and such suppres- 
sion effect is similar to previous studies about the impurity 
effect in the SDW state B25I1 . It is also notable that for a un- 
pinned vortex, the field-induced magnetic order is weak for 
the present parameters, and thus the suppression of the low- 
energy bound states is also expected to be quite weak. Since 
we used the same model and the same parameters, we actu- 
ally have a good consistent explanation for experimental re- 
sults in BaFe2_ x Co A As2 1231 12411 . Therefore, we conclude that 
the field-induced SDW cannot account for the disappearance 
of the bound states in Co-doped materials. 

We now study the low energy electronic structure for the 
pinned and unpinned vortex. The intensity plot of the LDOS 
at the real space with one pinned vortex and one unpinned 
vortex is shown in Fig. 3. As is seen, near the unpinned vor- 
tex the low-energy LDOS is enhanced and largest at the core 
center, indicating the existence of the low-energy bound state 
in the vortex core, agreeing qualitatively with the numerical 
results for hole-doped materials [6]. While for the pinned vor- 
tex with a strong impurity locating at the vortex center, the 
LDOS at the impurity site is zero. Near the center the spectra 
is suppressed significantly, and it seems that no bound state 
exists from the whole region of the pinned vortex. 

The LDOS spectra from the bulk to the vortex center are 
plotted in Fig. 4. The spectra without impurity are shown in 
Fig. 4(a). As is seen, there is a strong in-gap peak at the neg- 
ative energy for the spectrum at the vortex core center. The 
peak will spilt with a stronger peak at a negative energy when 
away from the center. The peak intensity decreases when fur- 
ther away from the vortex and finally the LDOS involves to the 
bulk one. All of the features are qualitatively the same with 
previous results for hole-doped samples [6], indicating that the 
changing chemical potential are not important to the low en- 
ergy electronic structure. The LDOS for the unpinned vortex 
with the impurity at the site (20, 10) is plotted in Fig. 4(b). The 
spectra are similar to those in clean systems, while the inten- 
sities of the resonance peak decrease. And some features out- 
side the SC gap appear, which are not a concern of the present 
work. The LDOS spectra for a pinned vortex with a impurity 
at (10, 24) are shown in Fig.4(c). For this potential the LDOS 
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FIG. 4: (Color online) The LDOS spectra for from the vortex center 
[the site (10, 14)] to the bulk respectively. Panel (a) is for the clean 
system and panels (b) and (c) are for the cases of one impurity with 
i m the impurity location. 



at the impurity site is zero (not shown here). Near the vortex 
core center, the LDOS is also suppressed signicantly and no 
low energy bound state exists for all the spectra we presented. 
This result is consistent with the STM experimental result re- 
ported in Co-doped materials 0. Actually, when away from 
the core center and when the SC coherent peak shows up, the 
coherent peak is enhanced at the negative energy. This is due 
to the impurity effect and discussed previously for the uniform 
SC state H] and SDW state JH] based on the same model. 
Since the enhancement is not significant enough and is near 
the gap edge, it may be covered up by the SC coherent peak. 
As a result, it may be difficult to be revealed by experiments. 

The LDOS spectra with different impurity potentials are 
shown in Fig. 5. For all the potentials considered, the low 
energy states are suppressed significantly. For the positive 
potentials, some out-gap features show up as the potential 
decreases, as seen in Figs. 5(a) and 5(b). For the moderate 
strength negative potentials (V s = -5), the bound states in- 
duced by the impurity is clearer at the negative gap edge. 
While such bound states reduce for stronger potentials ( V s = 
-8). The spectra for V s = -100 (not shown here) are almost 
the same as the spectra V s = +100. 
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FIG. 5: (Color online) Similar to Fig.4 but with different impurity 
potentials and i m = (10, 14). 



We here have given a reasonable explanation for the STM 
experiments on BaFe2- A Co A As2. Since the experiment is for 
the sample 5 = 0.1, the impurity concentration is much larger 



than the vortex density. For a isolated vortex and consider- 
ing the critical distance is insensitive to the direction, we can 
roughly estimate the pinned ratio as (1 - 0.9 1 "" 1 )- The ratio is 
99.5% for r c — 4 and 95% for r c — 3. Thus it is understand- 
able that the bound state cannot be observed in experiments. 

In summary, based on a two-orbital model and taking into 
account the impurity effect, we have studied the vortex states 
for iron-based superconducting materials . The vortex will 
be pinned as the impurity is close to the vortex core. The 
low-energy bound states show up for a clean system and 
could be suppressed by the impurity. Therefore, we have pre- 
sented a consistent picture for recent STM experiments on the 
Bai- A; K x Fe2As2 and BaFe2-. v Co A As2 and elaborated that the 
presence of the impurity accounts for the absence of the bound 
states in the vortex core for the BaFe2- A Co A As2 compound. 
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